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This paper presents a modeling and control of aerostructure developed by lattice-based 
cellular materials/components. The proposed aerostructure concept leverages a building 
block strategy for lattice-based components which provide great adaptability to varying 
flight scenarios, the needs of which are essential for in-flight wing shaping control. A 
decentralized structural control design is proposed that utilizes discrete-time lumped mass 
transfer matrix method (DT-LM-TMM). The objective is to develop an effective reduced 
order model through DT-LM-TMM that can be used to design a decentralized controller 
for the structural control of a wing. The proposed approach developed in this paper shows 
that, as far as the performance of overall structural system is concerned, the reduced 
order model can be as effective as the full order model in designing an optimal stabilizing 
controller. 


I. Introduction 

Earlier research studies conducted by NASA had shown that active control of wing twist and deflection 
at local wing sections can improve overall aerodynamic efficiency^® Two recent NASA projects to realize 
in- fight aeroelastic wing shape control are the Variable Camber-Continuous Trailing Edge Flap (VCCTEF) 
and Distributed Electric Propulsion (DEP) concepts. The VCCTEF strategy divides a trailing edge into 
many flap sections, each of which can be individually controlled subject to the constraints of an interstitial 
elastomeric interpolating surface. Chord- wise and span- wise adaptive wing shaping for various flight scenarios 
can be achieved in this manner . 4 DEP proposes to accomplish similar shape morphing modes, but with 
differential control of spatially distributed electric propulsion systems. 

The cellular composite materials proposed in this paper allow for high-performance structures to be 
incrementally assembled, modified, and repaired without requiring custom tooling. Fiber composites have 
been widely used to construct truss cores 5 and structural frames P These have all entailed the development 
of processes to continuously wind the fibers through a structure, and subsequently cure a resin matrix 
around them. Instead, the approach taken here uses parts with carbon-fiber beams and loops that are 
reversibly linked. These materials are assembled from mass-produced carbon fiber reinforced polymer parts 
with integrated interlocking connections. 

Modeling of high dimensional lattice structures using conventional FEM approach can be a challenge, 
mainly because it is difficult to analyze and visualize the integrated lattice structure in real-time. In this 
paper, we propose to use the concept of discrete-time lumped mass transfer matrix method (DT-LM-TMM) 
to model and analyze large dynamical systems. Due to the unique construction of the wing, as shown in 
Figure [l] we believe that the lumped mass approach will be appropriate with masses being located at the ribs 
and the spring being calibrated to the stiffness of the wingbox in between. The damping that is included is a 
combination of the extremely low internal damping of the carbon fiber and the friction due to the reversible 
attachment of the parts. The basic concept behind this approach was inspired by the work of Tan et al.P 
where the notion of modified transfer matrix method (M-TMM) approach was developed by utilizing the 
dynamic stiffness matrix of finite element. The primary goal there was to reduce the computational efforts 
involved in structural analysis. In this paper, in addition to incorporate the notion of M-TMM, we also 
utilize the numerical integration approach proposed in Kumar and Sankai® and develop a reduced order 
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Figure 1. Lattice-based composite cellular wing structure and its lumped mass model. 


discrete-time state-space model that is best suited for control synthesis. An optimal decentralized LQR 
controller is then designed for the reduced order model, and it is demonstrated that the overall structural 
system performance is comparable with that of the optimal LQR controller with full order model. 

This paper is organized as follows. In Section II we provide a brief overview of lattice-based cellular 
materials/components. Section III contains the introduction and some basics of discrete-time transfer matrix 
method, and the reduced-order model by following the proposed transfer matrix method is given in Section 
IV. Section V focuses on developing and formulating the decentralized controls, and Section VI contains 
numerical simulations and discussions. Concluding remarks are made in Section VII. 

II. Lattice-based Cellular Structures 

Cellular materials can incorporate multiple functions; these have conventionally been based on uniform 
periodic or stochastic structures.^ The discrete construction of digital composites introduces an independent 
design degree of freedom that allows global functional properties to be determined by the local placement 
of heterogeneous components. The lattice-based cellular structures will be assembled in a regular fashion 
from many discrete components, for which there is a small set of different types. These can be made and 
measured separately, and then combined in a hierarchical fashion, based on their discretization. Preliminary 
tests indicate that the macroscopic properties of an assembled structure can be accurately predicted from 
the behavior of its components, by using the topological connectivity of the lattice to define an equivalent 
mesh for modeling, and by using measurement of individual components to define the behavior of the mesh 
elements. The discrete construction of these materials will allow for combinatorial optimization of system 
designs. 
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Figure 2. Example of a lumped mass configuration for discrete-time transfer matrix method. 


III. Discrete Time Transfer Matrix Method 

For completeness we will provide a quick overview of the discrete time transfer matrix method provided 
by Kumar and Sanka]® before going into the adaptations for use in decentralized control. For a lumped mass 
model; see Figure [2j the equation of motion for a given subsystem n can be described by 

m n (ti)X n {ti) = t£(U) - T^(tj) + f n (ti ) (1) 

where r, ^ and r ^ denote the forces from right and left side of n, and f n the control force at n. Furthermore, 
can be described by 


x n = k n{U) [x%(ti) ~ X*(ti)] + C n {U) [x^{U) ~ (U) ,] (2) 

and (fc n ,c n ) denote the pair of stiffness and damping components at the left side of n. Because of repeti- 
tiveness of the construction of subsystems, the left side forces r % at n th subsystem are equal to the right 
side forces r ? f_ 1 at (n — l) th . Adjacently we know that there is no change in displacement across the mass. 
These conditions are represented as follows, 


T n = T n - 1 (3a) 

x n = x n (3b) 

The generalized form of acceleration and velocity can be represented as 

x n {ti) = A n (ti)x n (ti ) + B n (ti) (4a) 

x n (ti) = D n (ti)x n (ti ) + E n (ti ) (4b) 

where the format of (A n , B n , D n . E n ) depends on the chosen numerical integration scheme. Now, substituting 
§ and @ into § yields 

m n [A n x n + B n ] = T* - + fn (5) 

where 


T n = k n [ x n ~ x n- 1 ] + c n [( D n X n + E n ) L - (£>„_ 1 X n _ 1 + E%_ i)] . 

Furthermore, by combining and ( |3b| ), we can rewrite © in the matrix representation 
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which can be described in compact form as 

Vn = P nVn ■ 


(6) 

(7) 

(8) 
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The same process can be applied to (3a) and (J6| which results in 
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and therefore can be equivalently written as 


Vn = F nVn - 1 • 


(9) 


( 10 ) 


Equations and ( flQ| ) can be combined to propagate the ’’state vector” v spatially forward from the left 
boundary condition, and can be represented by 
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or equivalently as 
Similarly, Equation (3a) together with (J6| results in 
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or can be equivalently described as 


v n = H nVn+l ■ 


( 12 ) 


(13) 


(14) 


Figure [2] shows a sequence of lumped masses which are interconnected by stiffness and damping components 
modeled by spring and dashpot. The propagation of neighboring forces from both sides of the subsystem n 
is highlighted. 


IV. Modeling and Control of Lattice Structures 



Figure 3. The combination of left and right propagation. 


In this section, we make sure of the discrete-time lumped mass transfer matrix method to develop a 
control-centric reduced order model for control system synthesis. Figure [I] shows a section of aircraft wing 
built by utilizing the lattice-based cellular structure concept. The composition of lattice-based construction 
renders itself naturally to lumped mass system setup. Each section of airfoil provides a logical location of a 
lumped mass, and the connecting components between each airfoil can be modeled as spring and damper. 
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A. Relating node 0 to node n 

In this section, we work to develop a relationship between the nodal states of n and 0 (the first element on 
the left edge). We start from (3b) and apply it to § and ( flO| ) to render 


v n — PnF n V n _ i , 

and we can continue this process until we reach subsystem n, 

Vn = F n Q n Vo , 

where 


(15) 


(16) 


Qn = n PiFi 


(17) 


i = 0 


and Q n denotes the transfer function relating node 0 to node n. Note that Vq represents boundary conditions 
at the left edge. 


B. Relating node m to node n 


Last section we have propagated the nodal states from left to right and arrived at node n. 
we establish the relationship between node m and node n by propagating from right to left, 
from (3b) but this time we apply it to Equations (12) and (14) to attain 


In this section, 
Again, we start 


V n — H n +lJ n +l v n +l 5 


(18) 


and this process can be continued and eventually results in 

y R _ rp R 
u n ~ ± n u rn j 

where v ^ is the boundary conditions at right edge, and 


n 

Tn = n H i J i - 

i=m 


(19) 


( 20 ) 


and T n represents the transfer function from node m to node n. 


C. Combining left and right propagation 

The proposed DT-LM-TMM approach described earlier can be applied recursively, and by combining the 
process shown in Sections |A| and |B| yields 


Vn = FnQnvff 

v n = T nV * 

m(A n X n + B n ) = T* - + f n 


( 21 ) 


We can solve (21 ) for the state vector at subsystem n. Then, by treating as the right and left boundary 


conditions in Equations JT6| and (19), respectively, we can solve for each subsequent nodal states recursively, 


as illustrated in Figure [3[fdr the case of clamped-free boundary conditions. Given the state vector v the 


condition (3b) can be applied that allows the free-end condition on the left side, while the clamp is being 
applied to the right side. It should be noted that it is not necessary to solve for the internal forces if they 
are not of interest. 
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Figure 4. Five lumped mass system subject to a unit impulse. 


V. Decentralized Control Problem Formulation 


The DT-LM-TMM presented in Section [0 can be applied to generate full state dynamic structural 
responses, as was studied in Kumar and SankarP The simulation results will be presented in the next 
section. The focus of this section is to develop a control-centric reduced order dynamic model that is best 
suited for discrete-time decentralized control design. To this end, we only consider the lumped mass systems 
with clamped-free boundary conditions, which models an aircraft wing attached to a fuselage. We start with 


a small variation of the system of equations (21): 
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(22) 


m(A n x 2 + B n ) = T* - + f n 

By splitting the propagation to each side of the equations we are able to avoid some numerical issues 
and the problem becomes easier to solve. What is more important though is the explicit statement of B n 
and E n which we use to generate the previous states and formulate the discrete time decentralized control 


problem. In solving (22), we are able to attain an equation of motion that is based solely on the current and 


the previous conditions, as shown below 

C x rx ^ = CE n _ 1 E n - 1 + CE n E n 


+ G E ri j rl -E'n+l + Gr + Cl — TR2B n + f n 


(23) 
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Kumar and SankaJ^ provided a comprehensive list of numerical integration techniques, including Fox-Euler, 
Newmark (3, and Houbolt methods. In this paper, we make use of the Houbolt integration method because 
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its framework can be readily incorporated into formulating the discrete-time control problem. The Houbolt 
integration coefficients are given by 8 


An(ti) AT 2 

B n (U) = ~ [5aj(ti_i) - 4:x(U- 2 ) + x(ti- 3 )] 

D n (ti) = 6AT 

E n (ti) = -gXr[18a;(t,_i) - 9x(ti_ 2 ) + 2z(fi_ 3 ) 


(25) 


It is important to note that the Houbolt integration method is only valid after 3 time steps have passed. 
Kumar and Sankar provided an extended set of coefficients for the earlier time steps, but they are not 
particularly important from the controls point of view and thus ignored here. 


XlO' 
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Figure 5. Fast Fourier Transforms of the third mass. The circles and squares indicate the peak magnitude. 


Substituting (|25j) into ([23]) results in a discrete-time model which has smaller dimension, as shown below 
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+cr^(CE n _ 1 E n -i + CE n+1 E n+ 1 + Cr + Cl + f n ) 

x n 

which can be equivalently rewritten as 

X(U) = AX(t^ t ) + B(a + f n ), 

where the system matrices (A, B) and the exogenous input a are given by 
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(28) 


ot = C E n _ 1 E n —i + CE n E n + Cr + Cl 


(29) 
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Equation (27) is the equation of motion for the mass ra n , where a consists of the contribution of forces from 
neighboring masses and f n is the control input applied to mass m n . To design a stabilizing decentralized 
controller, we consider a to be a ’’weak” coupling between mass m n and its neighbors, and therefore we can 
treat mass m n as a local subsystem. 10 

It is important to note that T^, and G^i contain information on the configuration of the spring, mass, 
and damper systems that were present between subsystem n and the boundary conditions. The configuration 
of C x r and Ce u also change depending on the boundary conditions. This allows the decentralized subsystem 
m n to act as if it is directly connected to the boundary conditions through a series of combined mass, spring, 
damper systems. 


VI. Results and Discussion 



Figure 6. Comparison of total energy of the uncontrolled, continuous LQR controlled, and DT-LM-TMM 
decentralized controlled systems. 


To demonstrate the proposed DT-LM-TMM and the decentralized control concept developed in the 
previous section, here we consider the vibration control of a five lumped mass system with clamped-free 
boundary conditions subject to a unit impulse, as shown in Figure [4] The control force is applied at the 
third mass. 

Throughout this section, the parameters used for generating simulation results are given in Table [I] 
We first perform a frequency domain analysis for the lumped mass system using three different numerical 
integration schemes; namely, 3rd order Runge-Kutta, 4th/5th order ODE, and Houbolt. Figure [5] shows the 
Fast Fourier Transforms of the third mass subject to a unit impulse applied at the fifth mass for a duration of 
10 /is. The result indicates that, except some differences in magnitude, prediction for modal frequencies is the 
same regardless of the integration technique used. Furthermore, it is important to note that accumulation 
of numerical error using Houbolt integration method tends to be stabilizing, while two other integration 
schemes do not. Therefore, the Houbolt integration scheme was adapted into the proposed DT-LM-TMM 
approach. 

By applying the recursive substitution process from the left and right boundaries to the central third mass 
subsystem, we obtain a reduced order model of the form (27). For this reduced order model, we assume that 
the position and velocity of the third mass is available for feedback and design an optimal linear quadratic 
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Figure 7. The continuous LQR controller synchronizes the four right masses first and then regulates all masses 
to zero. 
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Figure 8. The DT-LM-TMM LQR controller regulates the controlled mass and to propagate to neighboring 
masses for suppression. 


regulator (LQR) feedback control. The weighting matrices Q and R are chosen as follows, 


Q = 


le 9 

0 

0 


R = 0.01. 


(30) 


It can be shown that the optimal controller developed for the reduced order system is able to suppress 
all the states. To highlight this observation, and for the purpose of comparison, we have also developed 


a continuous-time full state optimal LQR controller with the same choice of weighting matrices (30). The 


simulations for both full state and reduced order LQR controlled system are performed in MATLAB envi- 
ronment. Figure [6] shows the comparison of the total system energy for two approaches, and it indicates 
that they both are able to suppress overall vibrational motions with comparable performance. The average 
system energy over the simulation time for the full state continuous LQR is 8.923e _n and for DT-LM-TMM 
is 1.113e -11 . 

Though both LQR controllers are able to stabilize overall dynamics, the response of each mass behaves 
very differently for each feedback-controlled system. The full state continuous-time LQR controller tends 
to first kill off all the low frequency content and the primary modes by forcing the four right masses to 
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Table 1. Parameters for simulated lumped mass system. 


Parameter 
Number of Masses 
Mass 

Spring Stiffness 
Damping Constant 
Time Step 


Value 
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0.0323fcg 

1000000 ^ 

le -7 — 

ms 
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Figure 9. The control energy of the system for the full state continuous LQR and the DT-LM-TMM LQR. 


synchronize, then slowly regulates them to zero. This can be seen in Figure [7J On the other hand, the 
DT-LM-TMM LQR controller simply attempts to attenuate the position and velocity of the controlled mass 
(third mass) and this effect is then propagated through boundary elements to achieve overall suppression, 
as shown in Figure [8] 

The amount of energy that is needed to regulate the vibrational motions is an important consideration 
when comparing the performance of controllers. Figure [9] shows the control energy of the two controllers. 
The control energy is calculated by multiplying the force input and the distance the mass traveled in the time 
step that the force was applied. If the mass is moving in the same direction as the direction of applied force, 
then the input energy in Figure [9] is positive, it is negative if the mass is traveling in the opposite direction 
of force. We can see that with the full state knowledge the continuous LQR has many occasions where it 
is attempting to simply slow down the system (negative energy), while the DT-LM-TMM always utilizes 
positive energy. This suggests that DT-LM-TMM always utilize the control force f n along the motion of 
central mass or when it is greater than the exogenous input a. 

Figure [9] also shows one of the advantages of the DT-LM-TMM approach, in that we can see the full state 
continuous LQR controller tends to oscillate at very high frequencies, resulting in the nearly shaded in blue 
sections. In these simulations the highest frequency content of the full state continuous LQR is 1.37Mhz, 
whereas the DT-LM-TMM is about 1.6Khz. Though both are too high for any practical actuation purpose, 
but by selecting the time step of 10 /is we have set a hard cap on the bandwidth at lOKhz. We believe the 
frequency content, especially in the case of DT-LM-TMM, can be reduced substantially by simply loosing 
up the integration time step without compromising too much on the numerical accuracy. 

The analyses performed thus far were derived from the choice of weighting matrices given in (30). In 
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Figure 10. Closed-loop performance vs. control energy: (a) Pareto optimal curve for the reduced order model 
through DT-LM-TMM; (b) Pareto optimal curves for the whole system under the same set of simulations as 
in (a) and for the full state continuous LQR simulations. 


Figure [lOj we show the Pareto optimal curves of the LQR control design with the reduced order model 
through DT-LM-TMM over a range of Q and R selections. The plots offer an interesting and important 
control design trade-off, from which one may determine the “optimal” design point that renders the best 
closed-loop performance given the available control effort. We also compare the Pareto optimal curves of the 
DT-LM-TMM LQR controller with that of the full state continuous LQR controller and observe that within 
a reasonable range of Q and R both controllers render about the same level of performance. 

VII. Conclusions 

In this paper, we have confirmed the viability of the proposed DT-LM-TMM approach as a simulation 
technique, which provided a means of calculating the propagation of boundary conditions for lattice-based 
structures. Furthermore, the proposed approach renders itself naturally to be formulated as a decentralized 
control problem for a reduced order model of much lower dimension. A stabilizing controller utilizing LQR 
approach was then applied to this control-centric reduced-order model, and we observed that the performance 
was similar to, if not better than, the performance of a continuous-time full state feedback LQR controller 
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applied to the full order system. The advantages of proposed methodology were apparent, in that we only 
need to deal with a localized model of small number of states and need to feedback only local states, hence 
less sensor knowledge. In addition, we gain many advantages by utilizing the discrete-time approach, such 
as easy migration to flight control software and explicit control with maximum bandwidth. 
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